arXiv:1508.04239vl [cond-mat.str-el] 18 Aug 2015 


Typeset with jpsjS.cIs <ver.l.0> 


Full Paper 


Compact Representation for Specific Heat of Interacting Fermion Systems 
in Terms of Fully Renormalized Matsubara Green Function 

Kazumasa Miyake^ and Atsushi Tsuruta^ 

^Toyota Physical and Chemical Research Institute, Nagakute, Aichi 480-1192, Japan 
^Division of Materials Physics, Department of Materials Engineering Science, Graduate School of 
Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan 

(Received May 7, 2015 ) 

On the basis of the Luttinger-Ward formalism for the thermodynamic potential, the spe¬ 
cific heat of single-component interacting fermion systems with fixed chemical potential is 
compactly expressed in terms of the fully renormalized Matsubara Green function. 


1. Introduction 

It is technically important to obtain a closed-form expression for the specific heat of in¬ 
teracting or strongly correlated fermion systems in terms of the fully renormalized Matsubara 
Green function, which is inferred phenomenologically or is now feasible on the basis of nu¬ 
merical methods in one form or another, such as through dynamical mean-field theory^) and 
1/Ai-expansion theory.^) In the 1970s, studies on entropy and specific heat along this concept 
were published for interacting phonons^^ and normal Fermi liquids.Since then, a lot of stud¬ 
ies have been published in various fields not only in condensed matter physics®^ but also in 
nuclear physics®^ and quark matter physics.However, to the best of our knowledge, much 
still needs to be improved. For example, some contribution is missing in Ref. 3 as pointed 
out in Ref. 4, in which the relevant missing term was calculated from the perturbation in 
interaction. 

In this paper, for interacting fermion systems, we find that the relevant missing term 
can be compactly expressed in terms of the fully renormalized Matsubara Green function 
starting from the Luttinger-Ward functional for the thermodynamic potential. We show that 
the specific heat can be expressed compactly in terms of the fully renormalized Matsubara 
Green function. In Sect. 2, we explicitly derive an expression for the specific heat with 
fixed chemical potential /r for a single-component fermion system in the form 

X X “ P)Gka{ien) 

+ ^ E ^ \^ka{ien)GuMn)] , (1) 

(T k n 

where a stands for spin degrees of freedom, Np is the number of lattice points, £k is the 
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dispersion of fermions, and Gkaii^n) and are the fully renormalized Matsubara Green 

function and its self-energy, respectively. In Sect. 3, some remarks are given. 


2. Derivation of Formula for the Specific Heat 

In this section, we show how the specific heat is calculated if the Green function G{k, icn) 
of fermions is explicitly given. We discuss the case of the system whose Hamiltonian is given 
as 



where Nl is the number of lattice points. According to Luttinger and Ward,®’®^ the thermo¬ 
dynamic potential fl is given as 

fl = Ho + (3) 


where 


and 


cr n ^ k 


(4) 


H = - 


{In [1 - Gl^{ien)T.kaiien)\ +T.ka{ien)GkcT{ien)} , (5) 


where G^^{i€n) is the unperturbed Green function 


1 


( 6 ) 

^en - Sk + fj- 

and H' is the perturbation expansion for the thermodynamic potential in terms of the skeleton 
diagrams with the fully renormalized Green function G and the bare two-body interaction V. 
The explicit form of Feynman diagrams for H' is shown in Fig. 1, and its analytic form is 
given by 


cr,cr' '^2 h\h 2 

+ 2 E! ^ El ^ El El ^kicr{i(^ni)Gk2cr{'i(-n2)^CF(7{ki — ^ 2 ) 

a ni n2 L kik2 

1 ^ ^ 1 ^ ^ 1 ^ El El Gkia{i(-ni)Gk2a{i^n2)Gk2,a{'i-(^n3) 

a ni 712 ^3 ^ fci~fc3 


^Gfcl—fc2-|-A!3,(T ^€712 T *^n3)kjjtr(^1 ^2 )I^(t(^2 ^s) 

1 ^ ^ 1 ^ ^ 1 ^ ^~~1 ^ 1 Gk^cr{.'i(-n\)Gk 2 cr{i^n 2 )Gk 2 ,cr'{.'i(-n 3 ) 

a,a' ni 712 713 L ki^^ks 

^G ki — k 2 +k 3 ,cr'(j'^ni *£712 T *6713 ) [IGcr '(^1 ^2)] 

1 ^ ^ 1 ^ El ^ El ^ ^ 1 El Gkicri'>'^ni)Gk 2 ( 7 {i^n 2 )Gk 3 (T{i^n 3 )Gk 4 (T{i^n 4 ) 

a Til 712 ri 3 714 ^ fci~fc 4 
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0' = -iOAAAA^+i 

Fig. 1. Feynman diagram representation of the thermodynamic potential The part indicated 
by dots includes the remaining terms of the third-order perturbation and all the higher-order 
perturbation terms. 



xGk l—k3+k4,a(,'^^ni T 

xVaa{kl - k2)Vaa{k2 “ k‘i)Vaa{k^ “ ^ 4 ) +., (7) 

where the part indicated by dots includes the remaining terms of the third-order perturbation 
and all the higher-order perturbation terms. We should note here that the factor 2 in the 
denominator of each perturbation term of is not explicitly mentioned in the textbook of 
Ref. 9. Therefore, Vl'/T is given as 

cr,(j' Til n2 ^ k\k2 

^ ^ ^ \ _^ 

+ 2 'y ^ T"' y ^ y ^ y ^ Gk-^(j(ieri\)Gk 2 ai.^^n 2 )^aa{ki — ^2) 

(T ni n 2 ^ kik 2 

+ ^ Gk4aYni)Gk2aYn2)Gk3a'Yn3) 

a,a' ni n 2 ^ ki^k^ 

^G/-^_fc2+/c3,crK^^7i.i i^n2 ^^713) 

X {[V^a'ikl - k2)f - - k2)V,,,{k2 - fcs)} 

Tg y Y y ^^y ^^y ^^y 'y ^ Gkia{'i‘^ni)Gk2cr{i^n2)Gk3a{'i'^n3)Gk4a(j'^n4) 

a ni 412 ns n4 L k\r^k 4 

'^Gki—k2+k4,a{}^ni *£712 T '^^n4)Gk4—k3+k4,a{}^ni *£713 T ^£774) 

xVaa{kl - k2)Vaa{k2 “ k^)Vaa{k^ “ ^ 4 ) +.• (8) 

First, we differentiate this term with respect to an explicit temperature dependence other 
than that included in the fully renormalized Green function Gk^ {i(-nj) '■ 

= -^EEE E G'fc4<x(i£777)e*-i'Gfc,.(ie„,)e--V..,(0) 

a,a' 711 712 ^ kik2 

+ 2 E! E! E! E! ^fcicr (^£711)^7^20-(*£712 )Ki-(t (^1 ~ ^2) 

(T 711 712 ^ kik 2 

"7 2 ^ ^, ^ y y y ^ y ^ Gkiai'i^ni)Gk2aYn2)Gk3,a'i'^^n3) 

G^a' 711 712 713 ki^^k^ 

^ /S2 + fc3,cr^ (^^Tli i^n2 “t“ ^^713) 

X {[V^a'iki - k2)f - - k2)V^a'{k2 - fcs)} 


A 

&f 
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4e^^eeee E Gkia{i^ni )Gfi2crii£n2)Glf^cri'i'^n3)Gk^cr{i^n4j 


u—u 

ni 712 >T’3 ^4 ^ fci~fc4 


xGfc, 

—h2-\-h4^^a{i^n\ i^n2 “t“ ^^ria H“ ^^ 714 ) 

Xl4a(^l - ^2)V;7a(fc2 “ k^)Vaa{k3 “ ^ 4 ) +.• (9) 

^E^EE ^ka {ien)Gka{i(^n)- (10) 


2t2 


n k 


Note that all the differentiations in the present paper are performed with the chemical poten¬ 
tial // hxed. In deriving Eq. (10) from Eq. (9), we have used an explicit form of the perturbation 
expansion of the self-energy, as discussed in Refs. 8 and 9. On the other hand, with the use 
of the identity 


A 

&T 


An' 

T QT r2 


and Eq. (10), we arrive at the following relation: 


( 11 ) 


( 12 ) 


Since the second term of 12, in Eq. (3), is given by a compact form in terms of the fully 
renormalized Green function G and the self-energy S, the remaining problem in calculating 
the specific heat in terms of G and S is how to extract information included in fl'. According to 
the theory of Luttinger and Ward,^’®) the variation of 12' with respect to the fully renormalized 
Green function satisfies the following variational relation: 


,512' = ^ r J] A ^ Y.^^{ie^) 5 Gka{ie. 
a n ^ k 

Therefore, the entropy S' arising from 12' is expressed as 

, _ _A 

dr 


(13) 


A^ A TVi 


dGka{ien) on' 
dr dT ’ 


(14) 


where d/dr implies the total differentiation through the explicit T dependence of 12' and 
through those of G and S, while djdT implies a differentiation only through the explicit T 
dependence. With the use of Eq. (12), Eq. (14) is reduced to 


S'=-E’’Ev7E5^‘-(«- 


2Vl 


^ dT ¥ 


iEE^E Efccr (^£n)Gfcf7 (^^7i) • (15) 


k a n k 

Hereafter, for concise presentation, we abbreviate the spin variables cr in G and E, considering 
the case of a paramagnetic state. 

Then, the specific heat G' arising from 12' is given as 
AS' 


G' = T 


dr 
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^ Nl^ dT \ dr 

n ^ k ^ 


+ jVi 


dG/f (iCfi) 


dr 


_ f 1 /d!^'\ d /n'\ 1 ^ d 

|r(dr ^ V^y ^ 2 dr fc(*en)]| 




Nl 

n ^ k 


A. 

dr 


^ki^^r, 


dGk {i^r 

d^ 


+ ‘^'^kii^n) 


dG/j {i^r: 


d~ rp^k{i^n)Gk{iSn) “1“ [Efc(ien)Gfc (*6^,)] 


+ 9E^E Mr E AT i'^l^i^^^Gki^n)]- 


Nl ^ dr 

n ^ k 


(16) 


Note that (dfl'/dr)* in the second equality of Eq. (16) is the r-derivative of Q' only through 
the r dependence of Gk{i(-n) while d{VL'/T)/dT is that through the explicit T dependence of 
[Vl'/T) as in Eqs. (9) and (10). In deriving the last equality of Eq. (16), we have used Eqs. 
(10) and (13). For the discussion below, we define the first and second terms of Eq. (16) as 
G'l and G 2 , respectively: 




Nl 


A. 

dr 


Efe (fCri 


dGkjier, 

dr 


+ ‘^'^kij'^n) 


dGk (^^n) 

dr 


d~^'^kii^n)Gk{iSn) “1“ [Efc(ie^)Gfc(ie^,)] 


and 


C'2 = ^E]^E^E^ [^U^en)Gk.iien)] 


(17) 


(18) 


The contribution to the specific heat from the thermodynamic potential fl, given by Eq. 
(5), is calculated as follows. With the use of the definition in Eq. (5) and the relation 1 — 
G^{i(-n)'^k{i(^n) = G^{i€n)G'j^^{i€n), the entropy S from fl is expressed as 


Nl 


-E-E^E 


Gliie 
+ 


k L 

(iSfc(ze. 


1 dG^{ien) „ , NdG;;, (ien) 
+ Gk[^€n) -- 

dG k{i^ri 


dr 


dr 


’Gkii^n) T ^ki'if^n') 


dr 


Therefore, the specific heat G from Q is expressed as 


<5 = ^^ 

dr 


^E^E^E 


k L 


GJ(> 






dr 


dr 


(19) 
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d'^kii^n) ^ /. \ I V /" • \dGk{'i'^r, 
H 77^ Gk\i€-nj + 


dT 


dr 




iVi ^ dT 

n k 


1 dG^(ien) (*en) 


+ 


G0(ie„) dT 
d^k (*£n) 


Gkii(^n) “1“ ^ki^^n) 


' dT 

dGfc (i6fi) 


( 20 ) 


dT ^ ^ dT 

According to the definition of the self-energy S, i.e., G~^ = {G^)~^ — S, and = 

i^n — Cfc + /^) the following relations hold: 

dGj^ (^€72) _ X677 d^ki^^^n) 


dT 


dT 


and 


dGfc(ien) _ rr’O/'iV 'll' 


dT T 

With the use of these relations, Eq. (20) is reduced to 


( 21 ) 

( 22 ) 




Nl 


k 


'i I py /■ N r)^ /■ \dGk{i(-n) 

'~7^G^.{ien) -\ —Gk{ien) + - 


+yTy—yT— 

^ N dT 

a n ^ k 


r<0 / ■ \ I py /. \ I V z'• \dGk{i€n) 

7^Gf^{ien) -|- -^Gkyi^n) + Efc(zen) 


.(23) 


The first term in the brace of Eq. (17) and the last term in the second bracket of Eq. (23) are 
of the same size and have opposite signs, so that they cancel each other. The second term in 
the brace of Eq. (17) and the last term in the first bracket of Eq. (23) are also of the same 
size and have opposite signs, so that they also cancel each other. Therefore, by adding Eqs. 
(16) and (23), G( -|- C is expressed as 

^i + ^ = EE]^E|-;^ [TGk{ien)^k{ien)] + 2 ie„ [Gk{ien) - Gl{ien)] 


+T^ 


dT 


Itr. 


{Gkiien) - Gliien)) 


(24) 


Because of the identity G — G^ = G^GS and the relation d(ie„/T)/dT = 0, Eq. (24) is reduced 
to 

C’i + (5 = EE]^E|-;^ [TGfc(fe„)Efc(ie7,)] + 2fe„G^(ie„)Gfc(ie„)Efc(ie„) 


+Eien^ [Gfc(fen)Gfc(fen)Efc(fe„)] \ . 


(25) 


With the use of Eq. (22), the last term in the brace of Eq. (25) is transformed to 
T[Gfc(*en)Gfc(ie„)Sfc(ie„)] 


= [Gl.{ien)f‘ Gk{ien)'^k{i(^n) + T ienG^(fe„)— [Gfe(ze„)Sfe(fe„)] . (26) 
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Therefore, Eq. (25) is reduced to 

+ ^ = [TGk{ien)^k{ien)] + ie„G^(ieOG'fc(ie„)Sfe(i 


ler, 


-\-ienTG^{ien)-^ [Gk{i^n)'^k{i^n)\ 

+ienG^^{ien) [l - ie„G^(ie„)] Gfc(ie„)Sfc(ie„)} . (27) 

With the use of the relations ienG‘l{ien) = 1 + (e^ — fi)G^{ien) and 'i£n[G^{ien)]‘^ = 
—r[dG^(ie„)/dT], which is eqnivalent to Eq. (22), the right-hand side of Eq. (27) is rear¬ 
ranged as 


^^ + ^ = EE]^El-;i^ [TGk{ien)^k{ien)] + ienGl{ien)^ [TGk{ien)^k{^en)] 


dT 


dT 


—^e. 


7 ^) [^fe(^^n)] Gk(i£n)^k(j'^n)'^ 


an k 


EE A7r E 1 dr 


+ [l + “ ^AG^{i(-n)^ [rGfc(ien)Efc(ie,i)] 


dr 

+T{ek - r) Gk{ien)^k{ien) 

E E E “ 7^)G'fc(ien)E [rGfc(ie„)Sfc(ze„)] 


dr 


"^dr ~ ■ TGk{ien)^k{ien) 


EE^E ±[T{ek-fi)Gl{ien)Gk{ien)^ki 


NL^dT 

an k 


l^r 


(28) 


The thermodynamic potential Qq for the free electrons is given by Eq. (4) and is calcnlated 
straightforwardly as:®^ 




Therefore, Gq = —T{d‘^Qo/dT‘^) is given by 


1 




ek- fL 


glcfc-Ml/r 


[e(^fc-M)/^ -p l] 


2 • 


Einally, the total specihc heat G is expressed as 

C = Co + C+G[ + C^, 

where C( -I- C is given by Eq. (28). Namely, 


1 


d 


c = co + EvzEdx 

(T h 


T'^{£k - n)Gliien)Gk{ien)^k{ien) 


+Go- 


(29) 


(30) 


(31) 


(32) 
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Because of the relation G^GT, = G — G^, the second term of Eq. (32) is given by 


d 

dr 


h-)Gfc(ien) 


dr 


1 




(33) 


The last term of Eq. (33) is easily transformed as 


dr 


1 


^ h')G'fc(*eri 


iVu 

dr 


(-) 


tanh 


EiE 


efc - 


2^ 

{<^h-G/T 


2T 


[e(^k-G/T _g i]' 


= -Co, 


(34) 


where we have used the expression for Gq in Eq. (30). 
Therefore, the specihc heat [Eq. (31)] is reduced to 


c - A 

dr 




+1^2) 


(35) 


where we explicitly represent the specific heat as G^ to indicate clearly that it is obtained 
under the condition of a fixed chemical potential fi, and is given by Eq. (18). Equation (35) 
is compact and a natural extension of that of Gq, Eq. (34) or (30), and is simply the expression 
in Eq. (1). However, it is much easier to use Eq. (32) for numerical calculations because the 
second term Eq. (32) more rapidly converges than Eq. (35) concerning the summation over 
the Matsubara frequencies e„’s. 


3. Concluding Remarks 


Starting from the Luttinger-Ward expression for the thermodynamic potential, we have 
derived a compact and exact formula for the specihc heat (with the hxed chemical poten¬ 
tial /i) in single-component interacting fermion systems in the form of Eq. (32) or (35). An 
advantage of using these expressions is that the specihc heat G^{T) can be calculated at any 
temperature T in strongly correlated fermion systems provided that the single-particle Green 
function is given by some means, e.g., dynamical mean-held theory^) and 1/A^-expansion the¬ 
ory of a la Nagoya,^^ 

Since the above derivation is exact, the specihc heat G^{T) corresponding to to the lead¬ 
ing order of T in the low-temperature limit must coincide with the expression given by Lut- 
tinger.^^) Verifying this equivalence is left for a future study. 

The expression for entropy is given by integrating the relation {dS/dT)^ = C^{T) /T, with 
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C'^(T) given by Eq. (35), with respect to the temperature T as follows: 

E - f^)GkAK-,T') 

. cr ^ h n 

E dT' E ^ E E T')Gk„{ie'^-,T ')], (36) 

^ -^0 cr k n 

where e'^ = (2n + l)7rT' , and explicit T dependences (other than that through the Matsub- 
ara frequencies) of the Green function and the self-energy are expressed as Gka-{i^n',T) and 
^ka{ien]T), respectively. 

Although we have developed the formalism for calculating the specific heat on the ba¬ 
sis of the single-component fermion system with an instantaneous two-body interaction, the 
Luttinger-Ward formalism for the thermodynamic potential is known to be generally applied 
to other systems such as the interacting Bose systems and the electron-phonon interacting 
systems.®^ 
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